Redshift-drift in axially symmetric quasi-spherical Szekeres models 
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Models of inhomogeneous universes constructed with exact solutions of Einstein's General Rel- 
ativity have been proposed in the literature with the aim of reproducing the cosmological data 
without any need for a dark energy component. Besides large scale inhomogeneity models spher- 
ically symmetric around the observer, Swiss-cheese models have also been studied. Among them, 
Swiss-cheeses where the inhomogeneous patches are modeled by different particular Szekeres so- 
lutions have been used for reproducing the apparent dimming of the type la supernovae (SNla). 
However, the problem of fitting such models to the SNla data is completely degenerate and we 
need other constraints to fully characterize them. One of the tests which is known to be able to 
discriminate between different cosmological models is the redshift-drift. This drift has already been 
calculated by different authors for Lemaitre-Tolman-Bondi (LTB) models. We compute it here for 
one particular axially symmetric quasi-spherical Szekeres (QSS) Swiss-cheese which has previously 
been shown to reproduce to a good accuracy the SNla data, and we compare the results to the drift 
in the ACDM model and in some LTB models that can be found in the literature. We show that it 
is a good discriminator between them. Then, we discuss our model's remaining degrees of freedom 
and propose a recipe to fully constrain them. 
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I. INTRODUCTION 

The last two decades have witnessed a phenomenal 
increase of the available cosmological data. Analyzed 
in the framework of an FLRW homogeneous cosmology, 
they have yielded the Concordance model where more 
than 95% of the Universe content is of unknown nature. 
Among these 95%, around 75% of the Universe energy 
density is ascribed to the influence of some rather exotic 
component called dark energy. But dark energy has never 
been directly observed neither in the Universe, nor in lab- 
oratories, and it has very exotic properties, namely, it is 
a kind of fluid with negative pressure or a cosmological 
constant with an amplitude too small to account for the 
vacuum energy in the current standard model of particle 
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physics. This has led some authors to studying whether 
the observations could not be given a more natural expla- 
nation in the framework of inhomogeneous models con- 
structed with exact solutions of Einstein's field equations 
without any dark energy component. 

The models most often found in the literature are 
roughly of two kinds: one patch large scale inhomoge- 
neous Lemaitre-Tolman-Bondi (LTB) models, spherically 
symmetric around the observer [l|-KL7| . and Swiss-cheese 
models in which the patches can be either spherically 
symmetric LTB holes [TH-H] or non spherical Szekeres 
patches [Hm. 

Since the LTB solutions are determined by two free 
functions of the radial coordinate plus a radial coordinate 
choice, two independent sets of data are necessary and 
sufficient to define them, e.g . angular distance and galaxy 
number counts [lil . l27H30j |. In the most recent works 
devoted to solve the dark energy problem with zero-A 
LTB models, such solutions have been constrained by two 
or more sets of data measured on our past light cone [3| , 
0-@, [Hill, EMI, [IMll . By construction, we are 
thus left with a degeneracy as regards the ACDM model, 
because in both cases, homogeneous and inhomogeneous 
cosmology, the same data are reproduced without the 
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possibility to discriminate between the models. However, 
it has been suggested that this degeneracy can only occur 
if the LTB model is not smooth at the center; otherwise, 
the models are distinguishable [HI, HH . 

The degeneracy is even worse, when a few independent 
data sets are used, for models constructed with Szekeres 
solutions which are determined by five free functions of 
the radial coordinate plus a radial coordinate choice in 
the most general case. Even if this number can be re- 
duced to three in the case of axial symmetry, as we dis- 
cuss in Sec. [Vj the mere supernova data used to constrain 
the models in [251 ] is far from being sufficient for fully de- 
termining them. 

This is the reason why tests using the redshift-drift 
have been recently proposed to deal with this issue. The 
redshift-drift is the temporal variation of the redshift of 
distant sources when the observation of the same source 
is done at observer's different proper times in an ex- 
panding universe. It allows one to make observations 
on the past light cones of an observer at different cosmo- 
logical times, therefore giving access to a slice of space- 
time. It is thus a good discriminator between different 
cosmological models able to reproduce the observational 
data on our current past light cone. This effect has 
first been considered by Sandage [36| . then by McVittie 
[37} ■ In Friedmann-Lemaitre-Robertson- Walker (FLRW) 
models, when the Universe expansion decelerates, all red- 
shifts decrease with time. In FLRW models where the ex- 
pansion is recently accelerating, like in the ACDM model, 
sources with redshifts < 2.5 exhibit a positive redshift- 
drift. In [HI this effect has been proposed to test the 
"Copernican Principle" and to close the reconstruction 
problem of an LTB model from the luminosity distance 
as inferred from the supernova data. Other authors have 
since examined this effect for one patch zero-A LTB mod- 
els [sllHij. It has been shown in [4l[ that this effect is the 
only currently proposed one able to test, in principle, the 
LTB one patch models against the ACDM model "outside 
the past light cone" , provided spherical symmetry is but a 
mathematical simplification and one considers LTB mod- 
els as exhibiting an energy density smoothed out over an- 
gles around us. However, other proposals designed to test 
specific LTB models considered as a single exact space- 
time and relying on conditions inside the observer's light 
cone can be found in the literature, e.g. tests using the 
BAO scale fl6l . l3~J] . the kinematic Sunyaev-Zel'dovich ef- 
fect 32. 33] or the Compton y-distortion of the CMB 
spectrum [Hill]. 

Here, we want to study the redshift-drift in Swiss- 
cheese models which seem to represent more closely our 
observed Universe with its voids and filamentary pat- 
terns. 

For voids, the spherically symmetric shape of the holes 
can be considered as a good approximation since non 
spherical voids evolve towards a spherical configuration. 
Actually, it has been shown, using the top-hat void 
model, that the smaller axis of an underdense ellipsoid 
stretches out faster than the longer one implying that 



voids become increasingly spherical as they evolve [43|. 
Conversely, overdense ellipsoids tend to form pancakes 
and filamentary structures which are what we observe in 
the Universe. This is the reason why we choose to study 
a Swiss-cheese with non spherically symmetric overdense 
patches. 

An inhomogeneous exact solution of General Relativ- 
ity with no symmetry at all (no Killing vector) is the 
Szekeres model [44| . Among its sub-classes, the quasi- 
spherical Szekeres (QSS) case is best suited for our pur- 
pose since it possesses all three FLRW classes of models 
as an homogeneous limit and can therefore be matched 
to any kind of FLRW background. However, it has been 
shown in [45|,[46[ that in a general Szekeres model, generic 
light rays do not have repeatable paths (RLPs), i.e. two 
rays sent from the same source at different times to the 
same observer do not proceed through the same succes- 
sion of intermediate matter particles, implying thus an 
angular drift of the source on the sky. This implies that, 
in such models, the light rays emitted from a source at 
different time coordinates are emitted in different direc- 
tions and therefore reach different loci at the border of 
each patch. Since we integrate the null geodesic and 
redshift-drift equations for a sequence of such different 
time coordinates for each source, i.e. for each redshift 
z 1 this makes impossible the choice of the locus of the 
matching between two patches and thus impairs the con- 
struction of Swiss-cheese models. There are only two 
Szekeres classes, besides FLRW models, where RLPs ex- 
ist, the most interesting for our purpose being the axially 
symmetric Szekeres models, in which the RLPs are the 
null geodesies intersecting every space of constant time 
on the axis of symmetry. 

Actually, such a class of models have been used in 
[25| (hereafter BC) to construct QSS Swiss-cheeses which 
can reproduce the supernova data. The model which 
we consider here is BC's model 5 which best reproduces 
these data and which is a model with overdense patches 
matched to a open FLRW background. In this model, the 
observer is located at the origin of the first patch and the 
sources, at redshifts ranging from z = to z ~ 2.5, can 
be anywhere on null geodesies axially directed towards 
the observer. Of course, this Swiss-cheese model must be 
still considered as a very simplified toy-model since axial 
symmetry is not a generic property of structures in the 
observed Universe. 

Anyhow, it is important to compute the redshift-drift 
in such configurations and to compare it to the results in 
the ACDM model and other LTB models. 

Inhomogeneous exact models have been sometimes 
criticized on the ground that they exhibit more degrees 
of freedom than FLRW models and are thus able to fit 
more easily the data. Another interesting question to be 
addressed is therefore: how many independent data sets 
do we need to fully reconstruct an axially symmetric QSS 
model? 

It has already been shown in [38j that the combina- 
tion of luminosity distance and redshift-drift data allows 
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one to fully constrain a spherically symmetric spacetime. 
This applies, in particular, to LTB models. It has also 
been shown in [29j | that, for any given isotropic observa- 
tions of the apparent luminosity l(z) and of the galaxy 
number count n(z), with any given source evolution func- 
tions L(z) and rh(z), a set of functions determining a 
zero-A LTB model can be found to make the LTB obser- 
vational relations fit the observations. 

We demonstrate here, as a first step towards a more 
general theorem to be applied to QSS models without 
any symmetry, that an axially symmetric zero-A Szekeres 
model needs three independent observation sets to be 
fully reconstructed. 

The structure of the present paper is as follows. In 
Secim we present the Szekeres models and the particular 
QSS subclass used in this paper. Section [TTT] is devoted 
to the derivation of the differential equation giving the 
rcdshift-drift. In Sec. IIV( we display our result for the 
computation of the drift in BC model 5 and we compare 
it to that of the ACDM model and other LTB models 
studied in the literature. In Sec. [V] we discuss the issue 
of closing the reconstruction of our type of model with a 
sufficient number of data sets, and we propose a solution 
for axially symmetric QSS models. In Sec. lVIl we present 
our conclusions. 



II. SZEKERES MODELS 

The Szekeres solutions [Jj are the most general solu- 
tions of Einstein's equations one can obtain with a dust 
gravitational source. They have no symmetry, i.e., no 
Killing vector, and are therefore well-suited to describe 
a lumpy universe. Their metric in comoving coordinates 
and synchronous time gauge is 



ds 2 



z 2 dt 2 



■2a 



dr 2 



,2£* 



(dx 2 +dy 2 ), 



(2.1) 



where a and f3 are functions of (<, r, x, y) to be determined 
by the field equations. 

There are two families of Szekeres solutions. The class 
II family, where /?' = (here the prime denotes deriva- 
tive with respect to r), is a simultaneous generalization of 
the Friedmann and Kantowski-Sachs models. Its spheri- 
cally symmetric limit is the Datt-Ruban solution whose 
physical interpretation is not clear [U |H| . The class I 
family of solutions is obtained when /3' ^ 0. They con- 
tain the LTB solution as a spherically symmetric limit. 
We choose therefore this class of solutions to study the 
redshift-drift in Szekeres models. After solving the Ein- 
stein equations, its metric can be written, after a change 
of parametrization more convenient for our purpose [49( , 



„2_^+2 (<f'-^E'/E) 2 _ dr2 q>2 (d X 2 + dy 2 ) 



ds 2 = c 2 dt 2 - 



- k 



E 2 



(2.2) 



where e = 0, ±1, $ is a function of t and r, k is a function 
of r, and 
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(2.3) 
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with S(r), P(r), Q(r), functions of r. 

A. Quasi-spherical Szekeres models 



As it appears from (|2.2p . only e = +1 allows the so- 
lution to have the three Friedmann limits (hyperbolic, 
flat and spherical). This is induced by the requirement 
of a Lorentzian signature for the metric. Since we are 
interested in the Friedmann limit of our model which 
we expect to become homogeneous at very large scales, 
i.e., that of the last-scattering, we focus only on the 
e = +1 case. It is called the quasi-spherical Szekeres 
(QSS) model. 

Its metric, obtained with e = +1 in eq. (l2.2l) . becomes 
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ds 2 = c 2 dt 2 



1 - k 



E 2 



where 
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(2.4) 



(2.5) 



Applying the Einstein equations to the metric (|2.4[) 
and assuming the energy momentum tensor is that of 
dust, the Einstein equations reduce to the following two: 



c 2 



2M 



fc + ^A$ 2 , 



(2.6) 



where the dot denotes derivation with respect to t, A 
is the cosmological constant and M(r) is an arbitrary 
function of r related to the density p via 



npc 



2M' - QME'/E 
$ 2 ($' -§E'/E) '■ 



where k = 8nG/c 4 . 
The 3D Ricci scalar is 



; k (<I>k'/k-2<f>E'/E 



1 



(2.7) 



(2.8) 



The Weyl curvature tensor decomposed into its electric 
and magnetic part is 



E a f3 = C a lPd u^u 



M(3$' - $M'/M) 
3$3($/ _ $E'/E) 



diag(0,2,-l,-l), 



Ha/3 = ^V-gta-y^C'* 1 ' PSU'U 5 = 0, 



(2.9) 



where e a7/JJ/ is the 4-dimensional Levi-Civita symbol. It 
can be easily noticed that (|2.6p - (|2.9|) reduces to the 
LTB equations once E'/E is set to vanish. 
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Since (|2.6I) is the same in this model as in the LTB 
model, the bang time function, ts(r), follows from (|2.6[) 
in the same manner as 



dt 1 $' - <f>E'/E 

dz ~ l + z& - $E'/E' 



(2.13) 



d$ 



-fc + 2M/$- 



= c[t-i B (r)]. (2.10) 



All the formulas given so far being covariant under co- 
ordinate transformations of the form f = g(r), this means 
that one of the functions k(r), S(r), P(r), Q(r), M(r) 
and tsir) can be fixed at our convenience by the choice 
of g. Hence, each Szekeres solution is fully determined 
by only five functions of r. In the following, we choose 
S, P, Q, M and ts, and we make the coordinate choice: 
${t ,r) = r. 



B. Axially symmetric QSS models. Null cone 
equations. 

We have seen in Sec. Q] that only axially symmetric QSS 
models possess RLPs and are therefore more suited for 
an easy study of the redshift-drift in Swiss-cheese mod- 
els. Moreover, in these models, the only RLPs are ax- 
ially directed null geodesies. Since, in our Swiss-cheese 
models, we use only radially directed light rays, this im- 
plies, as shown in [51] (see also 52]), that the Szekeres 
model should be axially symmetric and that, accordingly, 
the rays we consider for the computation of the drift arc 
RLPs. The simplest axially-symmctric Szekeres model 
obeys 



III. THE REDSHIFT-DRIFT EQUATION FOR 
AN AXIALLY SYMMETRIC QSS MODEL 

The redshift-drift is the redshift increase or decrease a 
comoving observer looking at the same comoving source 
on her past light cone can measure while her proper time 
elapses. That means that the redshift of the source is 
measured on the observer's two different past light cones. 



O'(ro,t + 5t ) 



(z + 6z)(r,t + St) 




FIG. 1: The redshift-drift Sz of a source, initially at a redshift 
z on the past light cone of an observer at O, as measured 
by the same observer at O' after an elapsed time 5to of the 
observer's proper time. 



P(r) = xq = const, Q(r) = y = const. 

In this case the dipole axis is along x = xq and y = yo 
(or in spherical coordinates along the directions $ — 
and $ = — 7r). 

For the axially directed geodesies (da; = dy — 0), we 
obtain from the null condition in (|2.4I) 



dt _ 1 & - <PE'/E 
dr c y/l — k 



(2.11) 



The plus sign is for r e < r Q and the minus sign for r e > 
r , with r e , the radial coordinate of the source and r D , 
the radial coordinate of the observer. Since we put the 
observer at the origin r a = 0, we use the minus sign in 

dun). 

The redshift relation in this case is 1521 



ln(l + z) : 
or equivalently, 

dr 



If $' - <f>E'/E 



(2.12) 



VT^k 



An equation giving the redshift-drift in LTB models 
was derived in [40j. We adapt here the reasoning to 
obtain such an equation for the axially symmetric QSS 
model. 

We consider an observer O located at the origin, with 
coordinates {to,r = 0). After an elapsed time Sto, this 
comoving observer is at O' (to + 5io, r = 0). A comoving 
source, with radial coordinate r has, on the past light 
cone issued from O, redshift z(r) and time coordinate 
t(r). The same comoving source has, on the light cone 
issued from O', redshift 



Z(r) = z(r) + 5z(r), 
and time coordinate 

T(r) =t(r) +St(r), 



(3.1) 



(3.2) 



with t(r = 0) = t , z(r = 0) = Z(r = 0) = 0, Sz(r = 
0) = and St(r = 0) = St . 

From the geodesic equations (|2.13[) . the equation for 
the redshift is 



dz 1 + z $E'/E 



dz _ 1 + z $' - $E'/E 
dr c y/l — k 



(3.3) 
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Eq. (|3.ip can be written 

dSz(r) dZ(r) 



dr 



dr 



dz(r) 
dr 



(3.4) 



Inserting (I3.3[) in (13.41) and considering 6z(r) and St(r) 
as negligible with respect to z and t, we keep only the 
first order terms in Sz and St and obtain 



dSz 
dr 



l + z 



<i>'(t,r) ~ if>(t,r)E'/E 



6'(t,r) - $(t,r)E'/E Sz 



St (3.5) 



Now, (|3.2p can be written 

d<ft(r) dT(r) dt(r) 



dr 



dr 



dr 



(3.6) 



Inserting (I2.11[) with the minus sign in (I3.6P and con- 
sidering (5t(r) as negligible with respect to t, therefore 
keeping only the first order term in St, we obtain 



dSt(r) _ I 4>'(i, r) — r)E'/E 
dr ~ c y/l-k 



St 



(3.7) 



We consider the case where the redshift z is monoton- 
ically increasing with r and we replace the independent 
variable r by z. By using 



_d_ 

dr 



dz d _l + z&- $E'/E d 
dr dz c 



c y/l - k dz 
in (|3.6p and after some rearrangement, we obtain 



dSz 
dz 



&-$E'/E" 'l + z 
Now, we insert (|3.8I) into (|3.7p . that gives 

St 



(3.8) 



(3.9) 



dSt 
~cL7 



l + z' 



(3.10) 



which we integrate from the observer O at (to, z = 0) to 
the source at (t, z) to obtain 



St = 



St 

l + z' 



We insert this expression for St into 
the equation for the redshift-drift, 



dz V 1 



Sz 



1 $' - $E'/E 
(1 + z) 2 $' - <&E'/E 



Sta 



(3.11) 
and obtain 

(3.12) 



Once Sz is obtained by numerically integrating (|3.12l) for 
a fixed value of Sto, the redshift-drift follows from its 
definition z = Sz/Stg. 



IV. COMPUTATION OF THE 
REDSHIFT-DRIFT IN BC'S SWISS-CHEESE 
MODEL 5 

A. The Swiss-cheese model 

As seen above, in the case of axially directed geodesies 
the equations which describe light propagation and 
redshift-drift simplify significantly. Moreover, density 
fluctuations (12.71) and curvature fluctuations, (|2.8p and 
(|2.9p . are the largest along the axial axis. When a light 
ray passes through the origin E'/E —> —E'/E [52l. |53T|. 
Since, in the above equations, E'/E is always multiplied 
by $ or <!> which are zero at the origin, there is no discon- 
tinuity here. Now, a much smaller discontinuity in E'/E 
(which probably appears due to the imperfect match- 
ing in the BC Swiss Cheese model, and would probably 
not occur for a model with perfect matching), of order 
4 x 10 -6 Mpc -1 (to be compared with the larger discon- 
tinuity of order 10~ 3 Mpc -1 at the origin) occurs also at 
the boundaries. However, its magnitude is not sufficient 
to impact noticeably the results for the redshift-drift as 
can be seen in Fig. [5] 

When constructing a Swiss-cheese model, one needs to 
satisfy the junction conditions at the patch borders. In 
our models the Szekeres inhomogeneities (the patches) 
are matched to the Friedmann background (the cheese). 
These Szekeres patches are arranged so that their bound- 
aries touch wherever a light ray exits its neighbor. Thus 
the ray immediately enters another patch and spends no 
time in the Friedmann background. The matching condi- 
tions across a comoving spherical surface, r = constant, 
are: the mass inside the junction surface in the Szekeres 
patch is equal to the mass that would be inside that sur- 
face in the homogeneous background (this mass is defined 
as the mass energy density p integrated over the spatial 
volume inside the junction surface at a fixed time coor- 
dinate t); the spatial curvature at the junction surface 
is the same in both the Szekeres and Friedmann regions, 
which implies k$z — kpr 2 and (ksz)' = Ikpr; the bang 
time and A must be continuous across the junction. In 
our model, A = in both the patches and the cheese. 

Besides matching the inhomogeneous patches, one also 
needs to take care of the null geodesies. However, since 
we only consider here axial geodesies, the junction is triv- 
ial and requires only matching the radial, or equivalently, 
the time component [54j . 

The Swiss-cheese model we consider exhibits patches 
that are described by the following functions: 



M = 



PbKC* 



dr r 2 



S p exp 



(4.1) 



where pf, = f2„, !^r, Sl r 



ms^G, "m - 0.3, H = 68 km s 1 Mpc \ 
180 Mpc, t B = P = Q = and 
S = (10 3 + f) ±0 - 8 , with I = r/kpc and where - is 



S p = 0.6 and a 
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for propagation from the origin with E' / E = S' / S = 
-0.8/(10 3 + I) and + towards the origin with E'/E = 
S'/S = 0.8/(10 3 + £). 

As can be shown from (f2~7j) - ([23j) . for r > 400 Mpc 
this model becomes almost homogeneous. We join in- 
homogeneous patches at r — 400Mpc. By construction, 
this is not a perfect matching. Only the functions M and 
k are continuous along this boundary but not for exam- 
ple the curvature [eqs. (12.8[) and (|2.9j) ]. In BC's paper, 
this non perfect matching is investigated with the use 
of a so-called minimum void scenario which exhibits the 
same arbitrary functions as model 5 but with only one 
inhomogeneous patch, and therefore no matching. The 
result is that the x 2 s for reproducing the SNIa data are 
exactly the same in both models, the Swiss-cheese model 
5 and minimal void model 5, and are equal to 269 for 307 
measurements. This shows that the non-perfect match- 
ing has no impact at all on SNIa data fitting. Our re- 
sult for the redshift-drift, displayed in Fig. [SJ shows that 
this non-perfect matching does not impair significantly 
the calculation of this effect (no visible discontinuities of 
the derivative at the boundaries there) . This ensures the 
validity of our model for the computation of the redshift- 
drift. 

We place the observer at the origin r = of the first 
patch and we consider the light emitted by the sources 
in the first patch and directed towards the observer. As 
we have seen above, in this case E'/E is positive. After 
calculating the redshift-drift for all the comoving sources 
with < r < 400 Mpc, we are led to the second patch and 
consider the sources located between the border of this 
patch and its own origin. The light traveling towards 
the observer moves away from this origin and we have 
therefore to change the sign of E'/E. The same occurs 
when we are led to the other side of the second patch 
origin. Light emitted from there travels again towards 
an origin and the sign of E'/E needs once again to be 
inverted, and so on. 



and 

M 

t - t B (r) = -j— fc)372 ( sinh ^ ~ ( 4 - 3 ) 

where rj(t, r) is the parameter. 

2. Using the initial condition $(to,r) = r in (|4.2I) and 
substituting rj = t]q at t = to, we get k(r) as 

_ = M(,)(cosh, -1) 
r 

3. We obtain 770 by solving numerically (I4.3|) for t = to- 
The numerical value of to for our model is 11.63 
Gyrs. 

4. Then we find t(r) on the past light cone by nu- 
merically solving the null condition equation p.llj) 
for incoming geodesies. The sign of E'/E is chosen 
as positive(negative) depending on whether a light 
ray is moving towards(away from) the origin in the 
patch in which it is. 

5. Substituting t(r) in (|4.3|) we obtain rj(r) on the past 
light cone issued from O, using which in (I4.2[) we 
calculate $(t(r),r) and its derivatives on the past 
light cone. 

6. Then we numerically solve the first of the two equa- 
tions (|2.13l) for the redshift z(t(r),r). 

7. After having found z, we find the redshift-drift at 
this z by numerically solving (|3. 121) . 

The various functions obtained by this procedure are 
described in the next sub-section. 



C. The results 



B. The algorithm 

In order to calculate the redshift-drift, we proceed in 
the following manner: 

1. To find k(r) we proceed as follows: 

After substituting t — to and ts(V) = into (|2.10|) 
we fix a value of r, which is the upper limit of the 
integral, between and 400 Mpc. We want such 
a k which satisfies (|2.10j) . We find that to have 
a real value of the integrand over the whole range 
of integration k < 2M/r is necessary. But in the 
range < k < 2M /r, we do not obtain any root of 
(|2.10|) . Therefore we choose the case k(r) < for 
which the parametric solution is the same as in the 
LTB model and is given by 

M 

$(r,t) = — -rr (cosh 77 - 1) (4.2) 



The results obtained by the numerical computation de- 
scribed above are plotted in Figs. [2] [5] The radial coordi- 
nate extends up to about 4,500 Mpc (some 5.6 patches) 
and the redshift, up to 3. 

Figure [2] depicts the function E'/E over these patches, 
while the functions M(r) and k(r) are presented in Figs. [3] 
and S] respectively. Figure [5] depicts the redshift-drift 
for the Szekeres model and compares it with that for 
the ACDM modeL and for three LTB models: Alnes et 
al.'s void model [3[ and the so-called constrained GBH 
(cGBH) void model Q, both studied in (39|, and Yoo's 
hump model studied in |40l |. As is well-known, the 
drift for the ACDM model is positive up to some redshift 
38, 39] - this is of course a distinguishing feature for an 
accelerating FLRW cosmology. 

The redshift-drift for our Szekeres Swiss-cheese model 
is negative while for the ACDM model it is positive until 
a redshift 2.5. However, we wish to stress here that the 
minus sign of the drift appearing in our particular model 
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cannot be considered as a general feature of Szekeres cos- 
mology. This should actually be studied by an analysis 
extended to a number of other QSS Swiss-cheese models. 
This will be the subject of future work. 



EVE(r) vs r 
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FIG. 4: The k(r) function, determining the model evolution 
type, as a function of the radial coordinate r. 



500 1000 1500 2000 2500 3000 3500 4000 4500 

r(Mpc) 



FIG. 2: The E' /E(r) function as a function of the radial 
coordinate r. 
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FIG. 3: The effective gravitational mass M(r) as a function 
of the radial coordinate r. 
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FIG. 5: The redshift-drift (5z/Sto) as a function of the red- 
shift 2 for the axially symmetric QSS Swiss-cheese model, the 
ACDM model, the cGBH LTB void model 0] (courtesy [3§]), 
the Alnes et al. LTB void model [j| (courtesy : [39( |1 and the 
Yoo LTB hump model [U (courtesy : Hj) 



As stressed before, the redshift-drift in our Szekeres 
model is also compared, in Fig. El with that in three 
LTB models: two LTB void models, the Alnes et al.-type 
[H model I and the cGBH mod el B studied in [39| , and a 
LTB hump model described in [5fJ| whose redshift drift is 
shown in [4fjj. In the void models, an inner under-density 
makes a smooth transition to an outer region with higher 
density. In the hump model the density profile exhibits 
a large over-density near the center. In all three mod- 
els the observer is centrally located. These models are 
constructed in order to reproduce to a very good accu- 
racy different cosmological data sets among which the 



SNIa data. They show the same general shape for the 
drift curve as our Szekeres model does, with the drift be- 
coming more negative with increasing z. Note however 
that the magnitude of the drift in our Szekeres model is 
higher by a factor of about two, at a given redshift, in 
comparison to the LTB models (at those redshifts where 
the LTB curves show a decline with increasing z). We 
want also to stress here that, even though the redshift 
drift of an off-center source observed at the symmetry 
center of any LTB void model is always negative, this 
is not obligatorily the case in a LTB hump model [40| . 
The feature exhibited by the redshift drift in the par- 
ticular hump model shown in Fig. [S] cannot be therefore 



8 



generalized to the whole class of such models. 



V. FULLY CONSTRAINING THE MODEL 
DEGREES OF FREEDOM 

Inhomogeneous exact models have been sometimes 
criticized on the ground that they exhibit more degrees 
of freedom than FLRW models and are thus able to fit 
more easily the data. In Sec. IIV1 we have studied a par- 
ticular axially symmetric QSS Swiss-cheese model able 
to reproduce the supernova data and we have computed 
the redshift-drift for this model. However, another inter- 
esting question to be addressed is: how many and which 
independent data sets do we need to be able to fully re- 
construct an axially symmetric Swiss-cheese QSS model 
from background observations? 

Following 51], it has been shown in [52| (Theorem 3.1) 
that a constant- (x, y) null geodesic exists only in an axi- 
ally symmetric Szekeres spacetime 1 . 

One can always find a coordinate system for x and y 
where, in such a spacetime, the functions P and Q obey 
the set of equations 



(x ,y ) = (P,Q)=^P' = Q' = 0. 



(5.1) 



Hence, in these coordinates, the model is fully deter- 
mined, not by six arbitrary functions of r, as in the most 
general Szekeres model, but only by four. We choose here 
M(r), fc(r), t B (r) and S{r) or equivalently E'/E(r). 

Up to here, we have only at our disposal, to con- 
strain our model, the supernova data used in BC and the 
redshift-drift we have computed in Sec. Mil to be used 
for such a purpose in the future. The question is: how 
many independent data sets do we need to fully constrain 
the four arbitrary functions of r of such a model assum- 
ing these data are measured on the axial null geodesies 
directed towards the observer? 

It has already been shown in (38j that the combination 
of luminosity distance and redshift-drift data allows one 
to fully reconstruct a spherically symmetric spacetime. 
This applies, in particular, to LTB models. It has also 
been shown in [29] that, for any given isotropic observa- 
tions of the apparent luminosity l(z) and of the galaxy 
number count n(z), with any given source evolution func- 
tions L(z) and m(z), a set of functions determining a 
zero-A LTB model can be found to make the LTB ob- 
servational relations fit the observations. We can reword 



1 We wish to stress here a mis-statement in this Theorem 3.1. It 
is claimed that there is only one null geodesic with constant x 
and y in the axially symmetric case. This is not true: the null 
geodesies that pass through the symmetry axis may have any 
radial direction, and may cross the symmetry axis at different 
instants of the cosmic time t. So, this is actually a 2-parameter 
family of null geodesies, one parameter defining the azimuthal 
direction and the other defining the instant of intersection with 
the axis. 



this last statement as: the two arbitrary functions of r 
of a zero-A LTB model can be fully determined by two 
independent sets of observations realized on the observer 
past light cone. 

A natural guess is therefore that an axially symmet- 
ric zero-A Swiss-cheese Szekeres model needs three inde- 
pendent observation sets to be fully reconstructed. We 
demonstrate here, using a reasoning inspired from [29j | 2 , 
that this is indeed the case. 

Note that, in this section, we work with units in which 
c = G = l. 



A. Coordinate choice 

We denote the past null cone issued from the observer 
at (t = to, r = 0) and composed of all the axially directed 
null geodesies of the Swiss-cheese by t = t(r). From 
(|2~TTj) . t = t satisfies, 

dt _ $'[t(r), r] - $[t(r), r]E'/E _ S 7 - $E'/E 

ch ~ VT~^k ~ y/T^k ' 

(5.2) 

a hat denoting a quantity evaluated on the observer's 
past light cone. 

We now use our freedom to choose the radial coordi- 
nate such that, on this past light cone, 



$' - $E'/E 



1. 



(5.3) 



The equation for the incoming axial null geodesic is then 
t{r) =to-r (5.4) 
With this coordinate choice, (|2.7j) becomes 

(5.5) 



^ 9 M' - 3ME'/E 
yl-fc 



and, on the past light cone, (|2.6|) must be written, with 
A = 0, 



I 2 2M , 
$ = — k. 

$ 



(5.6) 



This equation possesses solutions in terms of a parameter 
r](t, r) depending on the sign of k. When written on the 
past light cone, they become 



K(r) 



M{r) 



(5.7) 



2 We wish to correct here an error in eq.(31) of [29^ . At variance 
with what is claimed there, the luminosity distance is not the 
same as the diameter distance but is related to it by the reci- 
procity theorem [55ll expressed by our eq. 115.171 1. 
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• For k > 0, 

K(r) = k (j) — 1 — cosf/ £ = rj — sinf/; (5.8) 

• For k = 0, 

*f(r) = l ?=(l/2)f f=(l/6)5f; (5.9) 

• For k < 0, 

K{r) = —k </> = coshr/— 1 £ = sinh rj—rj. (5.10) 

These expressions tell us the type of evolution for 
each region, based on the measured data: for k > 0, 
the model expands away from an initial singularity and 
then re-collapses to a final singularity (elliptic evolution); 
with k < the model is either ever-expanding or ever- 
collapsing, depending on the initial conditions (hyper- 
bolic evolution); k = is the intermediate case (parabolic 
evolution) . 

B. Differential equations for the arbitrary functions 

The total derivative of on the past null cone is 

dt 



d$ 
dr 

into which wc substitute 



dr ' 

and (15.41) to obtain 



(5.11) 



From 



E' IE = - 



we have 



d$ ~ r- 
dr 



(5.12) 



$ = ±4 



/2M 



(5.13) 



Since we are in an expanding Universe, we choose the + 

sign for <I> in (|5.13[) . Thus, we substitute (I5.13[) with the 
+ sign into (|5.12[) and obtain 



E' IE = i 




(5.14) 



This gives us a first order differential equation for the ar- 
bitrary function E(r) which can be very easily integrated 
once $ and the other arbitrary functions M(r) and k(r) 
are determined. 

We wish to write now a first order differential equation 
for M(r). We substitute ([5~T4]) into (|53|) . which gives 



dM 3M 
dr $ 



d$ 



>2M 



-3- + W— - k - Vl-k 
dr V $ 



47r ( o$ 2 v / l - fc. 



(5.15) 



Note that, at variance with the differential equation for 
M obtained in the LTB case in 29], (|5.15l) involves not 
only the observable quantity p, but also the unknown ar- 
bitrary function k. This will not invalidate our reasoning 
but make it a little more complicated. 

We aim now at deriving a differential equation for k in- 
volving observable quantities, the arbitrary function M, 
$ and its total derivatives with respect to r. 

The luminosity distance, Dl, of a source satisfies, on 
the past light cone of the observer, 



D L {z) = 



1 L(z) 
4tt£(z) : 



(5.16) 



with L being the absolute luminosity, i.e. the luminosity 
in the rest frame of the source, and i(z), the measured 
bolometric flux, i.e. integrated over all frequencies, emit- 
ted by the source at redshift z. 

The area distance, Da, also known as the angular di- 
ameter distance, is related to the luminosity distance by 
the reciprocity theorem [55j . 



D A {z) 



Dl{z) 
iX + z) 



2 ' 



(5.17) 



Therefore, once we obtain a differential equation involv- 
ing Da, it is easy to transform it, through (|5.17p . into 
a differential equation involving Dl, which is a quantity 
measured on the observer's past light cone, as shown by 

It has been demonstrated in BC that, for axial 
geodesies, Da can be written 



^ = - l -R a ,k^ D/ 
ds z 2 



(5.18) 



where k a is the null tangent vector, dx a /ds, and s is 
the null affine parameter. On an axially directed null 
geodesic k x — k v — and (|5 . 18[) can thus be written 

d ~ L)A = -\ [R u (k t ) 2 + 2R tr k t k r + R rr (k r ) 2 ] D A . 



ds 2 
We use 



ds 2 

With E,-= E 



d 2 D A d 2 D A 2 

(K ) 



dr 2 
- E„ 



dD A dk r 
dr ds 



(5.19) 



(5.20) 



E, ry = on the axially di- 



rected null geodesies, the geodesic equation for r becomes 



*1 ,-tur 



ds $i 



where 



$i k' 
$7 + 2(1 - k) 



(fc r ) = 0, (5.21) 



$! = $'- $E'/E. 



(5.22) 
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With our coordinate choice, on the light cone, k r = 
-fc*. Using (|5.2ip . we insert (|5.20p into (|5.19[) and obtain 



d 2 D A 
dr 2 



$1 



$! 2(1 -fc) 

SI = o. 



dD A 
dr 

(5.23) 



The Rtr Ricci tensor component vanishes here and the 
Rtt component is 



Rtt — 



$' - <f>E'/E 
and the R rr component can be written 



(5.24) 



1 



1 



{ , 2 ($ 2 - 1 + fc) (£'/£) 2 + 2^ 



- $'<t> 2 



<E>'(1 - fc) - fc' 



fc' 



(5.25) 



Substituting (|5.24[) and (|5.25p written on the light cone 
into (|5.23|) . we obtain a differential equation involving 
Da(t), dDA/dr and d 2 I?A/dr 2 . However, what the ob- 
server measures on her past light cone are not these quan- 
tities, but rather Dl(z) and its derivatives with respect 
to the redshift z. We apply therefore 



AD A dD A dz 



dr 



dz dr ' 



(5.26) 



and 



d 2 D A dD A d 2 z , d 2 D A fdz 

dr 



dr 2 dz dr 2 dz 2 
while also using (|3.3I) . written on the light cone as 



(5.27) 



dz $! 
dr VI - fc 



(5.28) 



and the L>l expression as given by (|5.1T|) . We obtain, 
after some calculations 



dD L 


2D L \ 


dz 


l + z) 


1 














D, 



4>i 



M 

VT - ! ^ 2(1 + z) vT - ! 

1) + + 2-^ - $E'/E 



l + z 



fc' n E' E.-? 2 , 
— - 2^^ $ + fc 



— ^ f 31^/^-2^-$' 



(*0 



= 0. 



(5.29) 



To obtain a second order differential equation for fc, 
we substitute into (|5.29[) E' /E, as given by (|5.14p . and 
(E'/E)' and the partial derivatives of $ as given in the 
Appendix. For readability purpose, we leave this task to 
the interested reader. These are functions of fc, fc', fc", 
M, p, dp/dz, z and the total derivatives of $ with respect 
to r up to second order. 

To be able to solve for M and fc we need a third differ- 
ential equation involving <I> and possibly M and fc. It is 
given by (|5 . 28[) into which we substitute the expressions 
for the partial derivatives of $ as given in the Appendix 
and that of E'/E as given by (|5.14j) . We obtain 



2M d$ 
$2 dr 




vT 



(5.30) 



We have thus obtained a set of three differential equa- 
tions for the unknown functions M , k and However, 
since the observables, p, dp/dz, Dl and derivatives with 
respect to z and z are given in terms of the redshift z, 
and not of the unobservable coordinate r, we need to 
write them as functions of r instead of z. This can be 
performed by using (|ST28|) and the solution z(r) of (fQ2|) 
given in below Sec. IV CI For readability purpose, we do 
not write down here the resulting equation and let the 
interested reader do the straightforward calculation for 
herself. 

We solve numerically the set of three coupled equa- 
tions (|5.15p . (15.291) and (15.30)) with the origin conditions 
M(r = 0) = k(r = 0) = 8(r = 0) = [|| and obtain 
M(r), fc(r) and $(r). Then, E(r) proceeds from (|5.14p . 

We have now at our disposal three arbitrary functions, 
M(r), fc(r) and E(r) which are sufficient for fully deter- 
mining the model with our coordinate choice (|5.3p . The 
last integration function t B {r) proceeds from M(r), fc(r) 
and $ as follows. 

We solve for rj from the first equation in (|5 . T[) and (|5.8[) 
- prT0]l , Then t B follows from 



t — t B — to 



t B 



M 



(5.31) 



C. Equation for z{r) 

Besides p and D^, another measurable quantity we 
have at our disposal is the redshift-drift as given by 
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(|3.12p . On the light cone, after replacing the partial 
derivatives of $ by their expressions given in the Ap- 
pendix and using (|5.28l) . this equation can be written 



fz_ 

dr 2 



1 



(1 + 5) 



dip 
dz 



dz 
dr 



0, (5.32) 



where zq(z) = (5z/5to)(z) is the redshift-drift for a given 
Sto as measured by the observer for a source at redshift 
z on her past light cone when she is at (to,f = 0). This 
equation gives a second order differential equation for z. 
With the origin condition z(r = 0) = by definition 
of the redshift, it can be numerically solved from the 
measurements of zq to yield z(r). This allows us to write 
the <E> total derivatives with respect to z instead of r, by 
using (f5T28|) . 



D. Algorithm 

To obtain the functions, M, k, E and ts from the 
galaxy number count, n(z), supernova luminosity dis- 
tance and redshift-drift observations, we propose to pro- 
ceed as follows: 

1. take the discrete data £(z,0,4>) and n(z,9,(f>), and 
correct them for known observational bias and se- 
lection effects. We do not consider here evolution 
effects. Average them over angles, fit them to some 
smooth analytic function of z and obtain i(z) and 
n(z). 



2. assuming a phenomenological L(z) law, use (|5.16p 
to obtain Dl(z) and then its derivatives up to sec- 
ond order. 

3. determine from the data a constant average galaxy 
mass m (we do not consider here galaxy mass evo- 
lution) and obtain p(z) — mn(z). 

4. assuming observations of the redshift-drift have 
been performed over an elapsed time Sto sufficient 
to provide a robust data set, z (z, 9, (f), average it 
over all angles, fit it to some smooth analytic func- 
tion of z and obtain zq(z) and its first derivative. 



5. determine z(r) by solving (15.32)) and use it to trans- 
form the observables z, p, dp/dz, Dl and its deriva- 
tives with respect to z into functions of r. 

6. determine M(r) and k(r) from the set of dif- 
ferential equations ([530]) . ([5715]) and (|5^9"|i . 

7. knowing M(r), k(r) and $(r), determine E(r) by 
integrating (|5.14|) . 

At this stage, we have obtained three arbitrary 
functions of r which, with our coordinate choice, 
are enough to fully determine the model. However, 
for completeness, we give below the recipe to com- 
pute the last integration function is(r). 



8. solve for fj from the first equation in (|5.7I) and (|5.8[) 

- 

9. determine i_e(r) from (15 .3 If) . 

This is the algorithm that could be applied to an ideal 
universe, where the redshift is monotonically increasing 
and where no shell-crossing is present to close the past 
null cone. We make such assumptions here but can be 
led to a more detailed study of the conditions imposed 
by them on the data in a future work. 

Note that this recipe implies also the existence of a 
number of particular directions of axial symmetry in the 
observable Universe. It constitutes another simplifying 
assumption made at this stage. 

To be complete, our above algorithm derivation should 
include a discussion of the existence and uniqueness of the 
differential equation set solutions. Given the difficulty of 
this task, we postpone it to future work and just assume 
here that solutions actually exist and that some of them, 
were they to be non unique, could be selected owing to 
their physical properties best designed to represent our 
Universe. 



E. Observables pertaining to the studied model 

The model studied in Sec. IIVI has already been shown 
in BC to be able to reproduce to a good accuracy the 
supernova data, i.e., the observed luminosity distance- 
redshift relation. We have also depicted in Fig [5] the 
redshift-drift for this model. To characterize it com- 
pletely, we need thus the mass density p on the past light 
cone of the observer at O, as shown in Sec. |Vj We have 
therefore computed it as a function of both the radial 
coordinate r (Fig. [5]) and the redshift z (Fig. [TJ. 

Notice that, while pj Pb possesses a Swiss-cheese-like 
feature on the spacelike hypersurface t = t (Fig. O, it 
becomes smoother on the past light cone, while exhibit- 
ing anyhow some cusps at the patch borders where the 
matching conditions are not perfect as already discussed 
in Sec. IIVI (Fig. [7]). We have been able however to fit it 
rather accurately to the smooth (1 + z) 5 5 curve (Fig. [7]) 
with maximum fractional error of 38% and mean frac- 
tional error of 20%. 

The redshift-drift appears smoother. Here, the non 
perfect matching is less visible. While all the model de- 
termining functions are constructed as in a Swiss-cheese, 
Fig. [5] giving the drift on the past light cone from O 
shows a very smooth curve. We have thus been able to 



fit it to the smooth az b 
a = -0.023,6 = 2.46, c = 



function with best fit for 



-0.029 and d = 1.2 with max- 
imum fractional error of 12% and mean fractional error 
of 2%. 

Hence, it is interesting to note that a very patchy un- 
derlying inhomogeneous model can correspond to smooth 
quantities measured on the observer's past light cone. 
Such a result might be an interesting outcome of algo- 
rithms of the kind proposed here, were the Universe such 
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p(t .*p b vs ' 




2000 2500 3000 
r(Mpc) 

FIG. 6: The ratio of the QSS model mass density p(to, r) and 
the background FLRW model mass density p;, at the initial 
instant to on the observer's world line. This quantity can be 
calculated but cannot be observed since it is in a spacelike 
relation to the observer. 



p(t ,z)/p b vs Redshilt 




FIG. 7: The ratio of the QSS model mass density p(to, z) 
and the background FLRW model mass density pt, on the 
past light cone of the observer at O, before smoothing (solid 
curve) and after smoothing as described in the text (dotted 
curve). 



that very inhomogeneous solutions able to represent it 
would correspond to smooth observables. 



VI. CONCLUSIONS 

The redshift-drift was recently proposed in the litera- 
ture as a way to discriminate between different universe 
models able to reproduce the current cosmological data 
without resort to any dark energy component [38l - l40| . 
It is considered as an interesting tool for getting rid of 



the degeneracy which can appear between different so- 
lutions. Its expression has been known for a long time 
in the Friedmannian framework and it was lately calcu- 
lated for the inhomogeneous spherically symmetric LTB 
models HJ. 

We have, in Sec. IIII1 generalized this calculation to 
the axially symmetric QSS model and given there a new 
analytic expression for the drift in this model. 

Then, we have applied this result to its numerical com- 
putation for the axially symmetric QSS Swiss-cheese BC 
model 5 which was shown in [25| to be able to repro- 
duce to a good accuracy the supernova data. Comparing 
it to the ACDM drift up to a redshift z = 3, we have 
shown that this effect can be a good discriminator be- 
tween both models since, (i) as it is well-known, the drift 
in the ACDM model is positive up to z = 2.5, while it is 
negative with an enhanced decline at high z in our model, 
(ii) at redshift z = 3, where the drift has become nega- 
tive for the ACDM model, its amplitude is much higher 
by a factor of ~ 14 in our model. 

We have also compared the drift obtained for our 
model with that displayed in [39|, H(| for LTB models. 
We have found that the redshift-drift in our model ex- 
hibits the same general shape as in these models, but 
that its magnitude is higher on average by a factor of 
about two at a given redshift. This has been discussed 
at length in Sec. IIVI 

However, inhomogeneous exact models have been 
sometimes criticized on the ground that they exhibit 
more degrees of freedom than FLRW models and are thus 
able to fit more easily the data. Therefore, in Sec. fVl we 
have addressed the issue of finding how many and which 
independent data sets we need to be able to fully recon- 
struct an axially symmetric QSS model from background 
observations and we have also derived an algorithm for 
implementing this goal. Of course, we can presume that 
the choice of the data sets able to constrain the model is 
not unique, but the observables we have used are the only 
ones which have allowed us to derive a detailed algorithm 
for dealing so far with the issue. 

Given the rather high amplitude of the drift in our 
model, some cosmologists might be tempted to use it 
carelessly as an actual discriminator between QSS Swiss- 
chccscs and other models. It was indeed studied in 
[HH how the drift could be measured in the future by the 
proposed E-ELT instrument CODEX or in [H, H3 b Y 
the VLT instrument ESPRESSO. These were claimed in 
[39l | to be able to discriminate between ACDM and differ- 
ent LTB models in the course of an observation decade. 
Therefore, since the drift in our model is higher by a fac- 
tor of two, some researchers might be lead to conclude 
that a five year observation period might be sufficient 
to complete a test for our model. Another possibility to 
discriminate between inhomogeneous and homogeneous 
Universe models by measuring the redshift-drift with the 
future space-borne gravitational wave interferometer DE- 
CIGO/BBO was put forward in [59|,|60( where the exper- 
iment duration was estimated to be around 5-10 years. 
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Now, we want to stress that the Swiss-cheese model 
we studied here is a mere toy model, not designed to be 
actually put to the test. We are not claiming that our 
Universe is constructed with axially symmetric patches, 
nor that we might be located at the origin of one of 
such patches or that the light emitted by distant sources 
should come to us following axially directed geodesies. 

However, we believe that cosmology is a very exciting 
science which is still in its infancy. For more than fifty 
years, homogeneous cosmology has been the main sub- 
ject of study and, despite a number of attractive achieve- 
ments, it is still facing huge unresolved problems, such 
as dark matter, dark energy, galaxy formation, etc. In- 
homogeneous cosmology is still younger and it needs to 
progress step by step. This is the reason why, even if 
after LTB models, Szekeres ones are now coming slowly 
into play, their tricky properties lead us to begin with the 
study of the simpler QSS models, hence axial symmetry 
and its mathematical simplifications. 

We are confident that the rather simple algorithm we 
display here to completely determine an axially symmet- 
ric QSS model from background observations will be gen- 
eralized in the future to the most general QSS model, 
allowing therefore to obtain a very interesting represen- 
tation of our Universe from cosmological data. We are 
determined to go on working in this direction. 



From (|5.11[) with our coordinate choice (|5.3[) and (fTJ) 

for $, it comes 



^7 d$ 2 
$' = — +$. 

dr 



(3) 



To obtain <£>', we differentiate (TTJ) with respect to r, 

then substitute $' as given by ((3]) and M' as given by 
(|5.15[) into the resulting equation, which gives 



where 



and 



2M 

$2 



2Md$ 



k' 



— + y/T^kB 

$2 dr 1 2 



2M , 
a — — k 



j3 = Airp® - 

$2 



(4) 



(5) 



(6) 



Now, we wish to calculate With our coordinate 

choice (|5.3p , we can write 
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$' = Vi~k + <$>E'/E, (7) 
which wc differentiate with respect to r to obtain 
k' 



$'E'/E + <P(E'/Ey. 



(8) 



The expression for <f>" is therefore obtained by substitut- 
ing into ([8j), the expression for $' as given by ([3|) and 
those for E'/E as given by (|5.14j) and for (E'/E)' as 
calculated below. 

Differentiating (|5.14j) with respect to r and rearranging 
yields 



Appendix 

We give here the expressions of the partial derivatives 
of $ needed to perform the calculations in Sec. [Vj 

The first derivative with respect to time is already 
known from (|5 . 1 3[) . We recall it here for completeness: 



$ = 



>2M 



-k. 



(1) 



To obtain $ on the light cone, we take the derivative 
of (|5.6|) with respect to the time coordinate which gives 



$ = -- 



M 

$2 ' 



(2) 



et\ _ i d 2 $ 

E J f d^" 
2M d$ 



1 

Ml 

$2 



1 d$ (1$ 

$2 dr \ dr 



$2 dr 



+ /3VT^k 



VT^k 



2 VSyr^fc V s 



After some calculations, we thus obtain for <!>'' 



d 2 $ 


1 




dr 2 


+ - 


3M 


1 


_ 2Md$ 




$ 


dr 



-d$ 

*dr~ 



k' ■. 



(9) 



(10) 
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Differentiating ([2]) with respect to r and substituting 
into the resuiting equation the expressions for $' and M' 
as given respectively by ([3]) and (|5.15[) . we obtain 



M ( d$ _\ , 8 

$3 \ dr / $ 



(11) 



Finally, to obtain <j>", we differentiate <&' with respect 
to r and obtain 



$" — $' following from our coordinate choice, with <fr" 

and $' as obtained from, respectively, (fTU|) and (j4} sub- 
stituted, and where we have also inserted M' as given by 
(|5.15|) and $' as given by ©. 



We see that this equation includes p' . We use (|5.28p 

and the expressions ((3]) of and ([3]) of <!>', to write p' 
in terms of dp/dz as 



p 2M d 2 $ 2M / d$ 
$ " = ~ TV + — 3 T" 

11M\ 1 



Vl - k ( 127rp- 



+ Vl-kB 



$3 y a 3 / 2 
2M d$ 



2M / d$ 
$2 I dr 

k' 




1 



$2 dr 2 
'3/3 d$ 



2M /d$\ , , 

— -7- +Vl~k[ + 4tt//$ 

<j>3 \dr / \ $ dr 



2M 



2M $ - 1 



a$ 2 



2y/T=k 2 



A-' 



— — h vl — kB — — 
$2 dr 2 



(12) 



where we have replaced d/dr(d$ /dr) by its expression 



P' 



1 + z 

'd$ 



2M d$ 



fc' 



$2 dr H 2 

fc I dp 



$ \ dr / $ \ dz 



(13) 



Examining the above expressions for (E' /E)' and the 
different partial derivatives of $ involved in the set of 
three coupled differential equations derived in Sec. [Vj it 
is easy to conclude that these equations relate $ and its 
total derivative with respect to r up to second order, the 
observables p(z) and its first derivative with respect to 
z, Dl and its derivatives with respect to z up to sec- 
ond order, and z, and the arbitrary functions with their 
derivatives, M, M', k, k', k", which allows us to solve 
for $(r), M(r) and k(r) with the algorithm describe in 
Sec. El 
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